Strong electron correlations in cobalt valence tautomers 
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We have examined cobalt based valence tautomer molecules such as Co(SQ)2(phen) using density 
functional theory (DFT) and variational configuration interaction (VCI) approaches based upon 
a model Hamiltonian. Our DFT results extend earlier work by finding a reduced total energy 
gap (order 0.6 eV) between high temperature and low temperature states when we fully relax 
the coordinates (relative to experimental ones). Furthermore we demonstrate that charge transfer 
picture based upon formal valence arguments succeeds qualitatively while failing quantitatively due 
to strong covalency between the Co 3d orbitals and ligand p orbitals. With the VCI approach, 
we argue that the high temperature, high spin phase is strongly mixed valent, with about 30% 
admixture of Co(III) into the predominantly Co(II) ground state. We confirm this mixed valence 
through a fit to the XANES spectra. Moreover, the strong electron correlations of the mixed valent 
phase provide an energy lowering of about 0.2-0.3 eV of the high temperature phase relative to the 
low temperature one. Finally, we use the domain model to account for the extraordinarily large 
entropy and enthalpy values associated with the transition. 



I. INTRODUCTION 

Transition metal complexes with redox active ligands 
have been the subject of extensive studies in recent years. 
Of particular interest are the valence tautomers based 
upon semiquinone and catecholate groups ligated to tran- 
sition metal ions . These molecules have intrinsic in- 
terest as candidates for mechanical and magnetic switch- 
ing devices activated by temperature, pressure, or irra- 
diation. They also present coordination environments 
and structural and magnetic conformation reminiscent, 
at least, of the local properties in certain allosteric met- 
alloproteins 

In this paper, we focus on the cobalt based valence tau- 
tomers of generic (high temperature) form Co(SQ)2Lp, 
where SQ represents a semiquinone ligand complex with 
formal charge -1 and spin 1/2, and Lp represents a neu- 
tral redox passive counterligand such as phen, py2X 
(X=0,Se,S,Te), tetramethylmethylenediamine (tmmda), 
or tetramethylethylenediamine (tmeda). When cooled in 
frozen organic solvents such as toluene, or in molecular 
solid form (typically under pressure), these molecules un- 
dergo a transition, in a narrow temperature regime (typ- 
ically between 100 and 400K), from a high spin, high 
volume form, to a low spin, low volume form [Q-^. The 
moment changes from order 4-6 /is at high T to 2 /is at 
low T, while the bond lengths of the inner coordination 
sphere typically decrease by about 10%, about 0.2A. 

Traditionally, it is assumed that there is a charge trans- 
fer from Co to the SQ ligands concomitant with the tran- 
sition, with the enhanced ligand field splitting arising 
from the bond contraction favoring low spin Co(III). This 
picture accounts for the observed susceptibility data, and 
is consistent with observed hyperfine splittings, optical 
absorption. 

However, there remain some puzzles in this accepted 
description of the valence tautomers, notably: 



• (1) On the theoretical side, density functional cal- 
culations PI have so far provided modest support 
to the above description, yielding stable high spin 
states for the high temperature geometry, and sta- 
ble low spin states for the low temperature geome- 
try. However, more Co electron charge is found at 
low T. 

• (2) For the Co(SQ)2(phen) complex, x-ray absorp- 
tion near edge spectra (XANES) spectra measured 
for the Co K-edge jl^, there is additional spec- 
tral weight at high energies (about 2-3 eV above 
the XANES peaks associated with transitions to 
3(i states) which has not yet been accounted for. 

• (3) The estimated entropy change AS* at the tran- 
sition inferred from fits to the effective moment is 
extraordinarily large, of order 10-15 R j5|, which is 
difficult to account for from softened vibrational 
modes at higher T. Concomitantly, a large en- 
thalpy of order 0.3 eV/molecule is required. Di- 
rect measurements of AS* for molecular solid 
forms give, on the other hand, a value of about 1.2 
R, consistent with the spin entropy change, which 
gives IS.H Ki 0.03 eV/molecule. 

Motivated by these experimental and theoretical puz- 
zles, we have set about to extend the description of the 
cobalt valence tautomers driven in large measure by a 
desire to assess the role of strong electron correlations 
in such systems. To this end, we have: (i) performed 
new, spin polarized, fully relaxed density functional the- 
ory calculations. These calculations provide an energy 
lowering of the high T-low T splitting which partly re- 
duces the disagreement between earlier DFT calculations 
and experiment. A Mulliken population analysis reveals 
0.2 fewer Co 3d electrons for the high-T form when com- 
pared to the low-T form, in marked contrast to the stan- 
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dard picture. We demonstrate, through an examination 
of the projected densities of states, that the arguments 
based upon formal oxidation states are quaUtatively cor- 
rect when one identifies the overall molecular orbitals de- 
rived from the nominal Cg symmetry about the Co site, 
(ii) We have developed a model Hamiltonian for the ac- 
tive electronic states which we solve by a variational wave 
function calculation which includes strong electron cor- 
relation effects on the Co site (outside DFT) . This calcu- 
lation amounts to configuration interaction with a physi- 
cally motivated basis set reduction which has proven use- 
ful for analyzing solid state systems. We find a significant 
further reduction in the high-T/low-T energy difference, 
with stronger electron correlation effects in the high-T 
phase (which has 0.3 less Co 3d electrons than the low-T 
phase). In particular, at high-T the high spin Co(II) sig- 
nificantly admixes with a high spin Co(III) state in which 
an SQ electron antiferromagnetically screens the Co mo- 
ment. We show that this a good quantitative account for 
the K-edge XANES data discussed above, (iii) We show 
that a domain model [^3|jlJ] can resolve the discrepancy 
between inferred and direct measurements of A5', 
Specifically, we show that existing experimental data is 
inconsistent with the assumption of a random mixture 
of high-T and low-T forms of the tautomer and is best 
explained by assuming the molecules form clusters of size 
~30-50. 

We stress that this description of the high-T phase of 
the valence tautomers makes them mixed valent in the 
sense used in the physics community to describe lan- 
thanide compounds such as CeSns or SmS, in which 
the ground state is a quantum superposition of states 
with predominantly two to three well defined valences for 
the rare earth ion. This mixed valence is distinct from 
the multi-site mixed valence of, e.g., the Creutz-Taube 
molecule. We distinguish this single site mixed valence 
from covalence, which is applicable to strong molecular 
orbital admixture at the single particle level. In this 
sense, a mixed valent single ion is intermediate between 
fully localized and fully covalent. Such molecular states 
have already been noted, e.g., in the case of cerocene 
where formal valence arguments generate the expectation 
of diamagnetic, tetravalent Ce ions, while observations 
yield temperature independent paramagnetism. The Cc 
ion is mixed valent in this case. 

We discuss the DFT results in Sec. II, the many body 
theory in Sec. Ill, the domain model in Sec. IV, and 
conclusions in Sec. V. 



II. ELECTRONIC STRUCTURE CALCULATIONS 

A. Computational Method 

We carry out the spin-polarized electronic structure 
calculations using the ab initio code SIESTA |16|,p7t. 



The Kohn-Sham equations are solved with the exchange- 
correlation calculated using the generalized-gradient- 
corrections(GGA) approximation, in the fully ab initio 
version of Perdew, Burke and Ernzerhof ||l8|. We use 
TrouUier-Martins norm-conserving pseudopotentials [Tg| ] 
in the Kleinman-Bylander form [ pO[ . For Co, we use non- 
linear partial-core corrections to account for exchange 
and correlation effects in the core region 111. The ba- 
sis set orbitals are obtained using the method of Sankey 
and Niklewski p^, g eneralized for multiple-!^ and polar- 
ization functions |23| . For the Co atom the basis set con- 
sisted of doublc-C 3d and 4s orbitals with 4p polarization 
orbitals. For C,N and O the basis set was double-C 2s 
and 2p orbitals whereas for H it was double-^ Is orbitals. 

The initial atomic coordinates are obtained using the 
experimentally determined high-temperature and low- 
temperature geometries for the phen complex. In both 
cases, following previous work |11|, the tertiary butyl 
groups were removed and replaced with H atoms (at a 
distance ~ lA from the C atoms). The supercell ap- 
proximation was used and a cubic supercell of dimen- 
sion 38 Bohr was taken to ensure sufficient vacuum be- 
tween neighboring tautomer molecules. The results ob- 
tained were checked for convergence by increasing su- 
percell size to 45 Bohr. Conjugate-gradient relaxation 
was performed on the system, with the experimental ge- 
ometries providing the set of initial coordinates, to de- 
termine the minimum energy configuration within DFT. 
The atoms were allowed to relax till the force on each 
atom was less than 0.04 eV/A. The electronic structure 
was determined both for the initial geometry and the fi- 
nal relaxed geometry and the results are presented in the 
following. 



B. Results 

The ground-state energetics obtained by us for the ex- 
perimental geometries are in agreement with previous 
ab initio calculations For the high-temperature ge- 
ometry, the electronic structure calculation yields high- 
spin (^ 3/2) Co coupled ferromagnetically to spin 1/2 
on each of the ligands. Furthermore we find that 
the low-spin configuration cannot be stabilized for the 
high-temperature geometry. In contrast, for the low- 
temperature geometry, the low-spin (~ 0) configuration 
is stabilized and the high-spin configuration is unstable. 
In the low temperature geometry there is a net spin 1/2 
for the complex which is located on the oxygen ligands. 
Furthermore, we find a net spin of 0.06 on the Co ion, 
consistent with the observed hyperfine splitting of 30G 
or less reduced by an order of magnitude compared 
to, e.g., metallic Co. Thus, our results are in agreement 
with experimental data which show a transition from a 
high-spin state to a low-spin state as the temperature is 
decreased. 
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The reason is of course, not complicated: the increased 
•id-2p hybridization upon contraction increases the hy- 
bridization induced hgand field splitting on the Co site, 
which disfavors Hund's rule alignment of the spins (the 
point charge contribution increases as well). The entropy 
cost is offset by the energy gain from increased ligand 
field energy. 

The ground-state energy difference between the two 
states without relaxation from the experimental geome- 
try is found to be 1.3 eV, which is in good agreement 
with the previously obtained value of 1.2 eV. However 
this energy difference is too high to be consistent with 
the room-temperature transition observed in the system. 
Experimentally the enthalpy change (A_ff) between the 
two states is inferred from fits to the magnetic moment 
to be ~ 0.3 eV/molecule § which is much lower than the 
theoretical value (as mentioned, the lone direct specific 
heat measurement suggests AH w 0.03 eV/molecule). 
However, even this experimentally inferred value for the 
enthalpy corresponds to an entropy change of ~ 14 i? 
which is unphysically large. Thus we are faced with a 
discrepancy between the experimental results and the en- 
ergetics obtained using DFT. 

This discrepancy can be partially resolved by not- 
ing that the experimentally determined ground-state 
geometry can differ slightly from the DFT ground- 
state geometry. Thus to allow the system to relax 
to the DFT ground-state configuration, we performed 
conjugate-gradient relaxation starting from the experi- 
mental geometries. Relaxing the geometries resulted in 
a substantial lowering of the energy difference between 
the high-spin and low-spin configurations to 0.6 eV, as 
compared to 1.3 eV for the unrelaxed geometries. In this 
relaxation process, the bond distances between Co and 
ligands changes very little, and thus the relaxed geom- 
etry might be argued as a better comparison point for 
experiment. In a previous study for the unrelaxed ex- 
perimental geometries, Adams et al j5| obtained a lower- 
ing of the ground-state for the high-spin configuration by 
performing a 'broken-symmetry' run. However, we have 
found that the broken-symmetry configuration is not sta- 
ble and upon relaxation the ground-state configuration is 
given by the ferromagnetic alignment of the ligand spins 
to the Co spin. However, the ground-state energy dif- 
ference obtained after relaxation is still too high to be 
consistent with the low transition temperature and it in- 
dicates that correlation effects which are not included in 
DFT play a significant role in determining the ground- 
state energies. 

Performing the MuUiken population analysis for the 
relaxed geometries leads to another discrepancy with the 
traditional interpretation of valence tautomerism. In this 
framework, the transition is described as a spin crossover 
phenomenon associated with charge transfer from the 
metal to the ligands. However the DFT results indi- 
cate that the net d-charge increases in going from the 



high-spin geometry (7.1) to the low-spin geometry (7.3). 
This should be contrasted with the traditional picture in 
which we should expect the net d-charge to be lower by 
1 in the low-spin geometry due to charge transfer to the 
ligands. 
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FIG. 1. Partial densities of states (PDOS) for Co-valence 
tautomers in the vicinity of the first ionization level, (a) Co 
3d-projected PDOS for the majority spin in the high temper- 
ature phase, (b) Co 3d- projected PDOS for the minority spin 
in the high temperature phase, (c) Co 3d-projected PDOS for 
majority spin in the high temperature phase. There is neg- 
ligible difference in this case between majority and minority 
spins, (d) PDOS jointly projected to SQ/CAT 2p and Co 3d 
levels in the low temperature phase. Note that the formal eg 
levels, for example, are strongly covalently admixed between 
the ligands and the Co ion. 

A resolution of this discrepancy is reached by analyz- 
ing the transition in terms of molecular orbital theory 
and studying the projected density of states (PDOS) for 
the two geometries. In terms of molecular orbital theory, 
the transition can be summarized as follows. The d or- 
bitals on the metal are split into the t2g and eg groups. 
The metal eg levels mix with the ligand a orbitals to form 
molecular orbitals which are bonding and antibonding in 
character. The bonding levels are expected to be ligand 
dominated and the antibonding levels metal dominated. 
The t2g levels don't mix with the a orbitals (ignoring, 
for simplicity, the mixing with the ligand tt orbitals) and 
the splitting between the t2g levels and the antibonding 
molecular orbitals (dominated by metal eg levels) is what 
corresponds to the crystal-field splitting. In the high tem- 
perature geometry the lower crystal-field splitting stabi- 
lizes the high-spin phase in which the antibonding eg or- 
bitals are occupied by two electrons and the t2g levels 
have 1 hole. In the low temperature geometry the Cg lev- 
els are unoccupied and the t2g levels are filled and the 
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remaining electron is transferred to the ligands. Thus we 
have a crossover from a high-spin configuration to a low- 
spin configuration associated with charge transfer from 
the metal to the ligands. 

A detailed analysis of the PDOS, shown in Fig. 1, for 
both the high-spin and low-spin geometries reveals that 
the above picture is indeed accurate, and yet it is also 
consistent with the DFT result that the net c?-charge is 
greater in the low T phase. This is explained as fol- 
lows : In the high T phase there is an occupancy of the 
eg — (7 antibonding molecular orbitals of 2 electrons in the 
spin-up channel. In the low T phase one of these elec- 
trons gets transferred to the Co t2g level and the other 
is transferred to the ligand SQ level. However this does 
not translate into a lowering of charge on the Co atom by 
one. There is significant covalency between the metal-Cg 
levels and the ligand-cr orbitals. So when 2 spin-up elec- 
trons are transferred from the Sg — cr antibonding orbitals 
that corresponds to a loss of just 1.4 for the Co atom and 
not 2. On the other hand because there is greater overlap 
in the low-T phase there is greater metal contribution to 
the occupied Cg — a bonding orbitals (which are ligand 
dominated). The difference between the metal contribu- 
tions to the Eg — cr bonding orbital is ^ 0.6 (higher in the 
low T phase). So the net metal eg occupancy in going 
from high T to low T changes by —1.4 -I- 0.6 = —0.8 i.e 
by less than 1 electron and not 2 as the formal picture 
indicates. This is more than compensated by the increase 
in the occupancy of the t2g level by 1 and hence we see 
that the net charge increases in the low T phase. 

Thus we have seen that the DFT results support the 
traditional qualitative picture for valence tautomerism at 
the same time clarifying that there is no effective charge 
transfer in the system. However there is still a significant 
discrepancy between the calculated and observed energy 
differences between the two states. In the next section 
we will explore the role of correlation effects in reducing 
this discrepancy. 



III. VARIATIONAL MANY BODY THEORY 

In this section, we show that electron-electron interac- 
tion effects can provide a further reduction of AH beyond 
DFT, and that they provide an explanation for hereto- 
fore unexplained spectral intensity in the high-T XANES 
spectra. Moreover, they confirm that like cerocene, the 
high-T phase of the valence tautomers is mixed valent in 
the sense used within the physics community to describe 
SmS, CeSna and other intermetallic compounds. 

A. Model Hamiltonian 

In order to investigate the possible role for correla- 
tion effects in the cobalt valence tautomers, we wish to 



consider the simplest model that possesses the salient 
features that are commonly accepted as essential for 
the combined mctal-ligand charge transfer (CT)/spin 
crossover of the tautomeric interconversion. We model 
this with an Anderson Impurity Hamiltonian p5| ] , in 
which local electron correlation effects are treated ac- 
curately on the Co site, but neglected on the 'metallic' 
carbon rings. 

We have taken the structure of the Co(phen) complex 
as the physical basis for our model. We represent the o- 
quinone ligands by their benzene ring skeletons and the 
oxygen atoms which are nearest neighbor to the metal 
site, and neglect octahedral symmetry breaking at the Co 
site. The model system therefore consists of the cobalt 
atom, the Sd-levels split by the cubic field arising from 
the local six-fold coordination of the oxygen atoms, the 
nitrogen atoms of the (N-N) complex, and the aromatic 
rings of the semiquinones. 

We restrict our model to the electronically active tt 
orbitals of the semiquinone ligands and the Sd t2g and Cg 
orbitals of the Co ion. 

The LCAO Hamiltonian using second-quantized nota- 
tion for the simple model we have just described is 



Ho 



(1) 



where c\^^ and Cjpfj are, respectively the fermion cre- 
ation and annihilation operators of the aforementioned 
orbital set where i and j are site indices, a and /3 are or- 
bital symmetry-adapted labels, and a refers to the spin 
degeneracy. The tiajp are on-site energies for kx = jf3 and 
are electron-hopping integrals for ia ^ jf3. The on-site 3d 
energies are 



where 



AeLF-- 



ed-ja — eda + AtLF 



{iaa I AV \ dja) 



eda 



ep{iacr) 



(2) 



(3) 



-0.4Ao if 7 



O.6A0 if 7 = xy, yz, xz 



where AV^ is the potential barrier that an electron hop- 
ping from the Co atom to one of the ligand N or O atoms 
must tunnel through, i is the site sum over the 4 oxygens 
and 2 nitrogen and a refers to the three 2p-orbitals on 
those sites. 7 sums over the degeneracy of irreps of Oh 
and eda- is the bare d-electron energy while ep(iaa) are 
the ligand atom on-site energies. Aq is the electrostatic 
contribution to the ligand field (LF). ed-ya- then plays the 
role of the on-site 3c?-orbital energy in the model. 

We include electron-electron interactions only for the 
Co 3d electrons. These are certainly the most significant 
owing to the more localized character of the states, and 
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experience with transition metal oxide solids show this 
to be a good starting assumption. Including also the 
spin-orbit coupling, we thus add to Tio 

tT>(T' 

+ Y.^m^-s, (4) 

i=l 

where Ud is the direct Coulomb integral and J/f < is 
the ferromagnetic Hunds' rule exchange coupling. The 
occupancy operator is defined as — d^^d^^ where 

d!y^ creates an electron in the Sd-orbital 7 with spin a. li 
is the orbital angular momentum of the ith d-electron and 
Si is the one-electron spin, ^(ri) is the free-ion parameter. 

For parameters, we make use of the tabulated on-site 
energies and also the distance parameterizations of the 
hopping integrals of Harrison . The hopping matrix 
elements are written as linear combinations of these in- 
tegrals using Slater-Koster theory |2^. Ud and Jh are 
treated as model parameters. We are currently working 
to constrain these parameters from information gained 
through ab-initio methods. However, we may obtain 
reasonable estimates by fitting the results of the model 
calculation to X-ray absorption data. This places con- 
straints on the charge transfer gap and also Udc-Ud, 
which is the difference between the Coulomb integrals 
describing the core hole-Sd electron and 3d-3d interac- 
tions, respectively. 

The high and low temperature expressions of the 
charge transfer gap are given by 

A^V~ ed - 0.4Ao - ^Jh + 6Ud - cl 
and (5) 
Aj^i^~ ed + O.6A0 + 6Ud -eL + 1.74eV 

where refers to the energy of an electron localized 
in a ligand orbital and 1.74 eV refers to the hybridiza- 
tion contribution to the ligand field splitting in the low- 
temperature geometry. 

Thus, the Hamiltonian is essentially that of a single- 
impurity Anderson Model, with the Co ion providing the 
localized states and the SQ ligands playing the role of the 
'metal'. Such a heuristic mapping to this model allows us 
to pursue a method of solution that we discuss in some 
detail in the next subsection. 

B. Variational Wave Function 

We proceed to solve for the ground state by perform- 
ing a variational configuration interaction (VCI) within a 
restricted basis set consisting of the lowest energy single- 
Slater determinants. This amounts to diagonalizing the 



Hamiltonian of Eqns.(l) and (4) within this basis of 
multi-electron wavefunctions. These determinants are 
built from the valence orbital set of our model. 

For the 'impurity' model described above, the varia- 
tional wavefunction method (VCI) has been shown to 
successfully describe transition metal or rare earth ions 
embedded both in metals and in metal oxides ||2|,|2^ . In 
the present molecular context we regard it as a physically 
motivated basis set reduction: we have 'divided' to con- 
quer by emphasizing first the strongest source of electron- 
electron correlations, expanding in the hybridization 
about the atomic limit for the Co ion. This said, the 
approach yields results which are non-perturbative in V, 
and systematically controlled by two handles: (i) the 
spin-|-orbital degeneracy (the lowest order results within 
a restricted Hilbert space become exact for large degen- 
eracy), with particle- hole excitations suppressed by in- 
verse powers of the degeneracy relative to the starting 
state, and (ii) the presence of excitation gaps which sup- 
press contributions from particle-hole excitations about 
the leading order configurations. 

The composition of the valence electron ground state 
will indicate the significance of correlations. The repul- 
sive nature of the Coulomb interaction will tend to drive 
d-electrons off the metal site. This should manifest as 
coherent quantum mechanical tunneling between states 
of Ud electrons and states belonging to the n^-l configu- 
ration (valence fluctuations). If such effects are insignif- 
icant, we should expect the type of single-determinant 
ground state that is treated very well within Hartree- 
Fock theory. 

The construction of the many-body basis is briefly out- 
lined. The states are Slater determinants that will consist 
of tensor products of cobalt states and ligand states, i.e. 
I rpSMs; La) where F is the crystal point-group irrep, 
/? is its degeneracy, S is the total spin, Ms is the spin 
degeneracy. La is the ligand state label which is used 
only if the a*'' orbital is occupied. Fig. 2 depicts the 
states we consider in the variational ansatz. We restrict 
our Hilbert space to the two low lying states containing 6 
or 7 Co 3d electrons; strong Coulomb repulsion legislates 
against other configurations. 

Magnetic susceptibility and x-ray absorption data, as 
well as first-principles calculations, have suggested the 
relevance of three closely- lying multiplets for the 3d- 
electrons: the ^Ti of the high-spin (h.s.) (5' = |), 
the 2e of low-spin (l.s.) (5" = i), and the ^Ai of 
l.s.-3(i® {S — 0). In the spin sector, we consider only the 
stretched state of maximal Ms within the S manifold. 
These would be expected to be the lowest-lying due to 
the Hund's rule energy. 

We also include an intermediate spin sector of states {S 
= 1, M^i^Ti and ^Tj) which couple to the h.s. and l.s. 
states through the spin-orbit interaction and to l.s.-3(i^ 
by CT. 

We assume that the SQ/CAT eigenstates of Hp are 
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adequately treated within the framework of Hiickel the- 
ory. This provides the standard approximate method of 
deaUng with tt* electrons. We employ these methods 
to find the unoccupied states of the two o-quinonc lig- 
ands and also their decomposition onto the constituent 
atomic orbital basis. Once these states have been ob- 
tained, Eqns.(l) and (4) may be diagonalized within the 
truncated basis of tensor product states. 
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FIG. 2. 



Singlc-Slatcr determinants 

4r 



which comprise the 
many-body basis. States (1-3) are '*Ti of h.s.-3ci'^, (4-6) be- 
long to ^T2 of h.s.-3d^ (7-12) are ^Ti, ^Ta of (S=l) 3d®, (13) 
is ^Ai which is Is.-SdP, and (14-15) is of l.s.-3d''. L repre- 
sents the doubly-degenerate spin-down LUMO of the quinone 
ligands. 



C. Results 

The Hamiltonian is diagonalized in both the high-T 
and low-T geometries using the optimized coordinates 
from the fully relaxed DFT runs. Fig. 3 is a histogram 
depicting how the quantum weight is distributed among 
the determinants of our variational ansatz ground state 
in both the high and low temperature phases. The high- 
T phase suggests a definite non-trivial role for correla- 
tions with a mixed valent state of 69% h.s.-Sd'', 28% 
h.s.-3d^ with an electron delocalized over the catechol 
ligand, and the remainder (~3%) in l.s.Scf . The single- 
determinant description works quite well in the low-T 
phase with 98.5% of the weight residing in the l.s.-3rf^ 
state, and (^1.5%) residing in the intermediate spin 3rf^ 
state. The small mixing between low-spin and higher- 
spin within a single configuration is due to the spin-orbit 
coupling. Therefore our low-T phase is exactly consistent 
with previous results of both ab initio calculations and 
experiment. 
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FIG. 3. Quantum weights of single determinants in 

Co(phcn) ground state wavcfunction. The vertical axis mea- 
sures \ ai \'^ where at are the variational coefficients. The 
model parameters correspond to = -0.05 eV and A^qj, 

= 0.42 eV. 

The most surprising aspect of these results is the reso- 
nance state in the high-T phase featuring large admixture 
with the h.s.-3rf^ state. There is a concomitant lowering 
of the energy such that the high-T phase is stabilized 
with respect to the low-T phase by ~0.4 cV. This en- 
tanglement of states would not be resolved by a mean- 
field theory like Hartree-Fock or DFT which both utilize 
a single-determinant groundstate composed of effective 
single-particle molecular orbitals. 

An estimate should be made of the contribution of 
these correlations to the high-T ground state energy. We 
can separate out the contribution to the resonant energy 
lowering due to covalency by finding the energy differ- 
ence between single determinant h.s.-3d^ and the high-T 
ground state of our calculation. This will be the energy 
omitted by DFT calculations. For completeness, we in- 
clude the possibility of the local electron tunneling to all 
excited states of the o-quinone ligands. This leads to 
an extra energy of ^0.2 eV that needs to be subtracted 
from the enthalpy difference between high- and low-T. 
This further reduces the 0.6 eV obtained from our re- 
laxed DFT calculations. The small magnitude of the 
gap Act seems to be indicative of the system's prox- 
imity to a charge-transfer instability within our model. 
This also justifies the high degree of inter-configurational 
admixture in the high-T phase. A^j, is nearly zero , 
slightly stabilized on the Ud = 7 side. A^^ is larger and 
has changed signs, so = 6 is stable. It will be seen 
in the next section that near-edge X-ray absorption data 
mandates some constraints as to where the high-T and 
low-T phases must reside relative to the instability. 
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D. XANES calculation 

We find that data taken in X-ray absorption studies 
near the cobalt K-edge corroborates the existence of 
the h.s.-3d^ /3cf mixed valence in the high-T state. We 
first discuss features of the experimental spectra. The 
on-set of the Is-Ap absorption edge is ~7712.2 eV. The 
region of immediate interest to us is the low-intensity 
pre-edge interpreted as being composed of quadrupolar- 
allowed ls-3d transitions. (Actually, the Co ions lack 
inversion symmetry, though if restricted to just nearest 
neighbors there is an approximate mirror symmetry for- 
bidding mixing between 4p and 3d states on site. The 
symmetry breaking by the SQ rings will yield an admix- 
ture of 4p and 3d which is likely to dominate the direct 
quadrupole matrix elements.) Data taken here provides 
a means to probe the nominal valence of the cobalt d- 
orbitals. As the temperature is reduced from Tamb, three 
main structures become apparent. The lowest energy fea- 
ture is a shoulder at 7708 eV that loses spectral weight as 
the temperature is reduced. In the high-T phase, there 
is an intense peak at 7709.6 eV (which shifts up in en- 
ergy by ^0.3 eV at lower temperatures) that gains ~56% 
in height. The third structure is a broad shoulder at 
^^7712. 3 eV that all but completely disappears at low 
temperatures. While the first two structures have been 
assigned the IsSd t2g and Is-Sd eg transitions, this latter 
feature has yet to be identified. 

We have used our model as the basis to calculate a 
theoretical XANES. It has been shown that the core 
electrons in transition metal and rare earth compounds 
do not efficiently screen the Is-core hole left behind af- 
ter the photoemission process. Thus the energies of 
the 3(i-valence orbitals are significantly renormalized by 
Coulomb interaction with the core-hole. Treated pertur- 
batively this amounts to a down-shift in energy for these 
orbitals. It has also been shown that in order to obtain 
good agreement with experimental data, a term must 
be added to the Hamiltonian that accounts for this core 
hole- valence electron interaction | |30[ |, 

5H = ^(1 - 4s^) d\,d^a (6) 

SO our Hamiltonian in Eqns.(l) and (4) becomes If 

H' = H + SH (7) 

where sj. is the Is-core electron creation operator, Udc 
is the ls-3c? Coulomb integral. The expression we use to 
calculate the X-ray absorption coefficient / for the ls-3c? 
transition is the following which is similar to a weighted 
many-body density of states 

lico) = ^ |(a' I T I vl/o) |2 SiLo - E^, + E,) (8) 



where lo is the incident photon energy, Ea' and | a') are 
the eigenstates of the Hamiltonian If which is the same 
system described by H in the presence of the core-hole 
perturbation. Eq and | ^'o) are the unperturbed ground 
state energy and wavefunction of our model calculation. 
T is the transition operator 

f = wY,d\„s^ (9) 

where we have neglected the orbital angular momentum 
dependence in W, the single-electron ls-3(i transition ma- 
trix element. 

We now discuss our results. Fig. 4(a) shows the spec- 
tra directly resulting from the model we have developed 
in the previous sections. The experimental XANES data 
has been superimposed over our results. Fig 4(b) depicts 
the result of a single-determinant l.s.-3d'' low-T ground 
state for the purposes of comparison; since the tautomers 
are on the edge of the charge transfer instability, it might 
be possible by suitable selection of the redox passive lig- 
and to access this low-spin Co (II) state. We have not es- 
timated the broadening generated by electronic coupling 
to vibrational degrees of freedom and other contributing 
factors to finite lifetime effects. This limits our compari- 
son to peak positions and relative weights. 

We first discuss the high-T result. Fig 4(a) features 
four main temperature dependent peaks that may be di- 
rectly related to the XAS data. We first note that our cal- 
culations reproduce the relative peak heights quite well. 
The peak at 7708.03 eV matches the shoulder found in 
that region and is the result of a down-spin Is electron 
transferring into the t2g. The shoulder feature has been 
previously assigned this transition. The peak at 7709.6 
eV is the consequence of the transfer of a down-spin elec- 
tron into an eg level, also consistent with previous as- 
signment. The energy distance between these two struc- 
tures ('--^l.BO eV) is approximately the value of the ligand- 
field splitting which we would expect on the basis of the 
assigned transitions. There are two very closely-spaced 
peaks at 7710.79 and 7710.8 eV. These are transitions of 
down-spin electrons into the feg shell of 3^*" final states. 
Going up in in energy by an amount equal to the ligand- 
field, we find a cluster of low-intensity peaks in the inter- 
val 7712.32 eV - 7712.42 eV that correspond to transfers 
into the eg levels of the 3ci'' final states. These latter fea- 
tures can be identified with the spectral weight centered 
about 7712.3 eV. This assignment is consistent with the 
disappearance of these peaks at low-T. To our knowl- 
edge, this is the first attempt to interpret these higher- 
energy satellite structures. This feature is thus unique to 
the type of highly admixed 3d''/3(i^ ground-state we are 
proposing for high-T. It can now be seen that the dis- 
tance between the peaks at ~7709.6 eV and the cluster 
centered about 7712.4 eV (and also the distance between 
the peak at 7708.03 and the peaks at ~7710.8 eV) of 
~2.5 eV provides a constraint for the value of Udc-Ud- 
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FIG. 4. (a) Model XANES calculation shown with experi- 
mental results for = -0.05 eV, A^y = 0.416 eV, E^mft = 
-1.05 eV, and Udc-Ud ~ 2.36 eV. Experimental curves show 
change in XANES with temperature. The solid line corre- 
sponds to the high-T phase and the ( — • • ■ — ) line is the 
low-T phase, (b) Low-T XANES calculation resulting from 
l.s.-Sd'^ groundstate. A'St = -0.004 eV, Eshift = 0.737 eV, 
and Udc-Ud = 2.43 eV 

We have introduced an energy shift between the high- 
and low-T spectra, Eshi ft ■ The pragmatic intent is to fit 
the data. The underlying physical reasoning is twofold. 
We have used identical model parameters for both the 
high-T and low-T calculations, allowing only the distance 
dependence of the hopping integrals to change. We would 
expect some renormalization of some of these parame- 
ters due to effects like screening. More significantly, we 
have omitted some terms from our simple model Hamilto- 
nian that would be quite different in the two geometries, 
e.g. cohesion energies, off-site Coulomb interactions, etc. 
Presumably, these would all have to be incorporated in 
such a shift. The low-T spectra in Fig 4(a) consists of 4 
main peak structures concentrated in the region 7709.66 
- 7710.8 eV. These are all associated with transferring a 
Is-electron into the empty Cg shell of the l.s.-3(i^ state. 
The heights and location of these agree quite well with 
the single broadened peak of experiment. The weights 
have not been altered relative to the high-T phase, giv- 
ing very good consistency with experiment. 

Fig 4(b) depicts a low-T XANES resulting from a l.s.- 
3(f groundstate. This regime should be investigated 
since this manifold of states lies in close energetic proxim- 
ity to the l.s.-3(i^ singlet state. As shown, this low-T state 
can be accessed with only very slight tuning of the pa- 
rameters. Such an interconversion would necessarily lack 



metal-ligand CT and may be relevant for some species of 
tautomer which only exhibits spin-crossover on the metal 
site. While the central peak at 7709.9 eV agrees well with 
the data, there is the feature of augmented weight at low 
energies (~7708.05 eV) at low-T rather than the suppres- 
sion seen in experiment. Since the absorption data seem 
to lack the twin peak structure (resulting from Hund's 
rule exchange splitting of up- and down- spin electron 
transitions to the singly-occupied eg shell), this seems to 
be an untenable possibility for the low-T groundstate, at 
least with current choices for redox passive ligands. 



IV. DOMAIN MODEL 

As noted in the introduction, another problem with 
the traditional interpretation of valence tautomerism is 
the extraordinarily large entropy change (10-15 R) asso- 
ciated with the transition. However, it should be noted 
that this entropy change is inferred only indirectly from 
measurements of the magnetic susceptibility and opti- 
cal absorption. On the the other hand, measurements 
of the specific heat as a function of temperature done by 
Abakumov et al jl^ , which constitute a more direct mea- 
surement of the entropy, give a much lower entropy value 
('--^ 1 R) . This discrepancy leads us to examine the set of 
assumptions used to infer the entropy change from the 
susceptibility data. These can be summarized as follows. 

At any given temperature T, the system is taken to 
be a two component mixture of the high-spin and low- 
spin forms of the molecule. The fraction of atoms in 
the high-spin form (/) is determined thermodynamically 
by minimizing the molar Gibbs free energy (G) which is 
given by 

G = fGhs + (1 - f)Gis - RT[f\n{f) + (1 - /)ln(l - /)] 

(10) 

where the last term is the entropy of mixing for a random 
two-component mixture. Minimizing the above equation 
with respect to / then gives us the functional form for 
the variation of / with T which used to infer the enthalpy 
(Ai7) and entropy (AS") change from the experimental 
data. 

However the above derivation makes the assumption 
that the system can be described as a random two- 
component mixture. If there is clustering present in the 
system, then the above equations have to be modified. 
The simplest generalization which accounts for this is 
the domain model |l^,|l^. In this model, clustering is 
accounted for by considering domains with n number of 
molecules per domain (n is like a mean-field parameter 
corresponding to the average size of the domain, with 
n = 1 corresponding to the random two-component mix- 
ture). In this case the entropy of mixing is now given 

by 
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5™. = -R/n [/ln(/) + (1 - /)ln(l - /)] (11) 
Correspondingly the variation of / with temperature is 
/ = l/[l + e^(i/T-i/T.)] (12) 

The key point is that the values inferred from exper- 
iment have used the above form with n — 1 giving un- 
physically large values for AH and A^. However if there 
is clustering (n > 1) then the inferred value of AH (and 
correspondingly A5') will be lowered by a factor of n 
leading to more reasonable values. This leads us to the 
question: How can we infer what the appropriate value of 
n is? According to the domain model the answer can be 
obtained from specific heat measurements which should 
indicate a peak in the specific heat (Cp) centered around 
the transition temperature Tc. The 'jump' in the specific 
heat is given by Cp{Tc) - i(Cz«(r,) + Chs{Tc)) where 
Cis{Tc){Chs{Tc)) is obtained by linear extrapolation at 
T = Tc of the low(high) temperature specific heat values. 
The number of molecules per domain n is then related to 
this jump by 

n = \§^\C,{T,) - liCisiT.) + CUTc))] (13) 

It should be noted that the specific heat measurements 
will provide a direct measurement of of AS (and thereby 
AH) and within the domain model this can be used to es- 
timate n. For the tautomers, only one set of experiments 
done so far have measured the specific heat. From the 
data given by Abakumov et al for the bpy complex, we 
make the following estimates : The 'jump' in Cp corre- 
sponds to ^ 150 J/mol K which gives a change in entropy 
of AS* ~ 1.2 i?. Using these values and the equation for n 
in the domain model we estimate n ~ 50. Thus the data 
for the specific heat shows evidence for clustering in the 
tautomers and is inconsistent with the original assump- 
tion of a random two-component mixture. The latter 
assumption leads to unphysically high inferred values of 
AH and AS* and we propose that measurements of the 
specific heat can be used to resolve this discrepancy and 
to self-consistently check for the applicability of the do- 
main model to the system. 

V. CONCLUSIONS 

In conclusion, we have extended the theory of cobalt 
valence tautomers to resolve some of the puzzles sur- 
rounding various data. In particular, we have shown 
that using fully relaxed density functional theory solu- 
tions goes halfway towards solving the mystery of large 
enthalpy values inferred from fits to optical absorption 
and magnetic susceptibility data. However, the relaxed 
DFT result of 0.6 eV is stiU off by a factor of 20 from 
the sole direct measurement of the enthalpy found in the 



tautomer literature |l^. We have also shown that DFT 
provides little evidence for net charge transfer off the Co 
ion to the SQ ligands (indeed, the 3d weight increases 
modestly in the low-T phase!). On the other hand, a 
careful examination of the PDOS shows that the qualita- 
tive picture emerging from the formal valence arguments 
is still correct: the orbitals with the largest admixture of 
3d — Cg symmetry are depopulated in going through the 
transition. 

Next, we showed that a variational configuration inter- 
action approach could produce a further lowering of the 
high-T enthalpy relative to the low-T entropy by produc- 
ing a 30% admixture of Co(II) and Co(III) configuration. 
At the lowest variational order, we estimate an 0.2-0.3 eV 
lowering of energy due to electron correlation effects be- 
yond DFT, which brings us within an order of magnitude 
of the directly observed enthalpy. This high-T state is in- 
trinsically "CI" : no unitary transformation in the single 
particle orbital space can remove the multiconfiguration 
character. It is, in fact, a mixed valent state associated 
with the Co ion and the more extended and weakly cor- 
related states of the SQ rings. We have also shown that 
this high-T state, when proper accounting is made for the 
core-hole interaction of the 3d electrons, can explain the 
previously observed spectral weight shift well above the 
dominant ls-3d peak for high temperatures. The weight 
is ascribed to transitions to the admixed high spin Co(III) 
configuration. 

Finally, we have proposed a domain model to account 
for the large AS*, AH values inferred from susceptibil- 
ity and optical absorption data (as compared to direct 
measurement). In this view, the transition corresponds 
to a broadened first order phase transition of molecular 
clusters of order 20-60. In the dilute solvent, it might 
be possible to observe such clustering with elastic light 
scattering experiments. Although the heat capacity mea- 
surements in the interesting temperature range are dif- 
ficult, we certainly would encourage further experiments 
on tautomers. 

We hope to apply this combined DFT and VCI ap- 
proach to other tautomer systems and to analyses of 
the closely related phenomena observed in the prosthetic 
complexes of metalloproteins and metalloenzymes. 
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